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Abstract 

The West Antarctic Peninsula (WAP) has been suffering an increase in its atmospheric temperature during the last 50 years, 
mainly associated with global warming. This increment of temperature trend associated with changes in sea-ice dynamics 
has an impact on organisms, affecting their phenology, physiology and distribution range. For instance, rapid demographic 
changes in Pygoscelis penguins have been reported over the last 50 years in WAP, resulting in population expansion of sub- 
Antarctic Gentoo penguin (P. papua) and retreat of Antarctic Adelie penguin (P. adeliae). Current global warming has been 
mainly associated with human activities; however these climate trends are framed in a historical context of climate changes, 
particularly during the Pleistocene, characterized by an alternation between glacial and interglacial periods. During the last 
maximal glacial (LGM~21,000 BP) the ice sheet cover reached its maximum extension on the West Antarctic Peninsula 
(WAP), causing local extinction of Antarctic taxa, migration to lower latitudes and/or survival in glacial refugia. We studied 
the HRVI of mtDNA and the nuclear intron fifibint7 of 150 individuals of the WAP to understand the demographic history 
and population structure of P. papua. We found high genetic diversity, reduced population genetic structure and a 
signature of population expansion estimated around 13,000 BP, much before the first paleocolony fossil records 
( — 1,100 BP). Our results suggest that the species may have survived in peri-Antarctic refugia such as South Georgia and 
North Sandwich islands and recolonized the Antarctic Peninsula and South Shetland Islands after the ice sheet retreat. 



Citation: Pena M. F, Poulin E, Dantas GPM, Gonzatez-Acuna D, Petry MV, et al. (2014) Have Historical Climate Changes Affected Gentoo Penguin (Pygoscelis papua) 
Populations in Antarctica? PLoS ONE 9(4): e95375. doi:10.1371/joumal.pone.0095375 

Editor: Gyaneshwer Chaubey, Estonian Biocentre and Tartu University, Estonia 

Received December 3, 2013; Accepted March 25, 2014; Published April 23, 2014 

Copyright: © 2014 Pena et al. This is an open-access article distributed under the terms of the Creative Commons Attribution License, which permits 
unrestricted use, distribution, and reproduction in any medium, provided the original author and source are credited. 

Funding: Funding was provided by: Instituto Ecologfa y Biodiversidad, Instituto Antartico Chileno (projects: M_04-10, T_27-10, G_06-1 1), Corporacion Nacional 
de Ciencia y Tecnologia (Master grant), Association of Polar Early Career Scientists, the National Institute of Science and Technology Antarctic Environmental 
Research (INCT-APA) (CNPq Process n : 57401 8/2008-5), FAPERJ (E-26/1 70.023/2008), FAPERGS (process 09/0574-7), WCS (processes 2008-5 and 2009-5), the 
Brazilian Ministries of Science, Technology and Innovation (MCTI), of Environmental (MMA), Inter-Ministry Commission for Sea Research (CIRM) and the University 
of Vale do Rio dos Sinos - Unisinos. The funders had no role in study design, data collection and analysis, decision to publish, or preparation of the manuscript. 

Competing Interests: The authors have declared that no competing interests exist. 

* E-mail: jvianna@uc.cl 



Introduction 

The West Antarctic Peninsula has been described as one of the 
areas most affected by global warming due to atmospheric 
temperature increment; it has experimented an increase of 5— 
6°C during the last 50 years [1—4]. The ecosystems of the WAP are 
principally regulated by ice dynamics, including its extension, 
periodicity and inter-annual variation [1]. Thus it is expected that 
changes in atmospheric and ocean temperatures will have a great 
impact on biological systems, affecting their phenology, physiology 
and distribution range [5,6]. At present there is strong evidence 
related to the effect of climate change on biological interactions 
and food webs [7] . Also, evidence of differential response between 
ice-obligated and non-ice-obligated marine mammals and birds 
has been found [5,8-10]. 

Penguins have been described as climate indicator species; 
changes in their populations over the last 50 years have been 
documented, mainly associated with changes in ice dynamics [1 1— 
1 5] . Among the Pygoscelis species, P. adeliae (Adelie) is the most ice- 
dependent species; its range reaches higher latitudes and has a 



circumpolar distribution. In contrast, P. antarcticus (Chinstrap) is 
found almost exclusively around the Antarctic Peninsula [16], 
while P. papua presents a distribution range including most of sub- 
Antarctic islands and the Antarctic Peninsula. Pygoscelis papua 
(Gentoo) populations have been increasing along the Antarctic 
Peninsula in the last 50 years; this has been mainly associated with 
an abrupt increment of temperatures in the Antarctic Peninsula 
that affects sea ice extension, which modulates access to breeding 
sites and the establishment of krill stocks. These changes in sea-ice 
dynamics have allowed southward colonization of P. papua 
populations [12,13,17-19]. Contrarily, populations of P. adeliae 
and P. antarcticus have been decreasing in the WAP, mainly related 
to the availability of krill [18,20]. 

Current global warming has been described by the Intergov- 
ernmental Panel on Climate Change (IPCC) without precedents 
over decades to millennia mainly associated with human activity. 
Nevertheless, this global trend is framed in a context of historical 
climate changes, particularly during the Pleistocene, characterized 
by an alternation between glacial and interglacial periods [21]. 
These warming and cooling events had an impact on the extension 
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and collapse of ice masses in Antarctic, periodically altering land, 
fresh water and marine habitats [22]. In the last glacial maximum 
(LGM—2 1,000 BP), Antarctic Peninsula ice sheet reached its 
maximum extension on the West Antarctic Peninsula (WAP). 
During the LGM decreasing temperatures and habitat loss were 
the main limitations for the permanence of biological communities 
[22]. The thickness of the ice sheet on the WAP, around 1.5- 
2.5 Km as well as ice shelf extension [23] caused a significant 
reduction of flora and fauna through local extinction of Antarctic 
taxa [24—26], migration to lower latitudes [27] and/or survival in 
glacial refugia [28]. There is also evidence that during past climate 
change that affected ice extension, penguin populations responded 
by abandoning their breeding colonies towards more profitable 
habitats [29], [30]. For example, as shown by paleoecological 
records, Antarctic P. adeliae penguins responded by abandoning 
nest colonies and dispersing to suitable habitats. Richie et al. [31] 
found the existence of two mitochondrial DNA lineages associated 
with different refuge histories during the LGM mainly attributed 
to the expansion of sea ice and loss of ice-free coastal areas 
required for nesting [29,30]. In contrast to Antarctic species such 
as P. adeliae, ice-independent species are expected to have been 
even more affected during the LGM, reducing their distribution to 
lower latitudes such as sub-Antarctic islands or to refuges near the 
Antarctic polar front. This has been supported by molecular and 
fossil evidence that suggests postglacial re-colonization of the 
Antarctic Peninsula by highly dispersive species such as marine 
mammals and birds [26,32], as well as shallow marine inverte- 
brates [33]. For example, molecular studies of Mirounga leonina 
show a post-LGM re-colonization of Victoria Land about 7,500- 
8,000 BP from Macquarie Island [32]. However, Emslie et al. [34] 
found that P. papua occupation of the Antarctic Peninsula did not 
occur until the late Holocene; the oldest record for this species is 
about — 1, 100 BP. Thus, we expected to find a clear recent signal 
of expansion for P. papua of less than 2,000 years. However, if the 
fossil found by Emslie et al. [34] was an outcome of the re- 
colonization after the LGM due to the expansion of one or more 
P. papua populations that were isolated during the Pleistocene 
climate change, we would expect to find an older signature of 
expansion and a marked genetic structure in the WAP region. 

In this study we used mitochondrial DNA (mtDNA) HVR-I and 
nuclear DNA (nDNA) $fibint7 sequences to evaluate past changes 
in the demographic history and population genetic structure of P. 
papua, to assess the dynamics of Antarctic populations during the 
LGM. Understanding how penguin populations have responded 
to past climate changes would allow predicting how future changes 
may affect their population dynamics. 

Methodology 

Sample Collection and DNA Extraction 

For penguins capture and handling, we used procedures with 
the least amount of stress. In order to minimize the disturbances 
we might be causing in the colonies, the penguins were captured 
several meters from the nesting sites when the animals left to the 
water. The penguins were captured using a hand-held net, and 
then immobilized manually as described by Wilson [35]. During 
the complete procedure the penguin remains with his head 
covered to reduce stress and the beak immobilized. Blood samples 
were collected from adult individuals from each studied colony. 
Samples were taken with 23 G needles from the brachial or 
external metatarsal vein (~0.5 mL) and stored in 96% ethanol for 
posterior analysis. After samples were collected, the individuals 
were released in the same place they were captured and they were 
monitored until finds a shelter or returns to the water. 



Pygoscelis papua blood samples were collected from five nesting 
areas in the Antarctic: Elephant Island (61°13'S; 55°21'W); 
Admiralty Bay (62°23'S, 58°67'W) and Ardley Island (62°13'S, 
58°54'W) on King George Island in the South Shetland Islands; 
and the Chilean Antarctic bases Gabriel Gonzalez Videla at 
Paradise Bay (64°49'S, 62°52'W) and Libertador Bernardo 
O'Higgins at Covadonga Harbor (63°19'S, 57°54'W; Figure 1) 
on the West Antarctic Peninsula. All samples were collected with 
permits in accordance to Annex II, Article 3 of the Protocol on 
Environmental Protection to the Antarctic Treaty, and the 
regulation from the Scientific Committee on Antarctic Research 
(SCAR) provided by the Chilean Antarctic Authorities (permits 
INACH 44/2012), and the Brazilian Antarctic Authorities 
through the PROANTAR and Environmental Ministry (N°21/ 
2011 and N°20/2012). The INACH and PROANTAR permit 
included authorization for sample collection for all locations and to 
develop research on Ardley Island and Admiralty Bay, the only 
locations for which a specific permit was required. Bioethics 
permit was provided by Universidad de Conception and Pontificia 
Universidad Catolica de Chile which was required to obtain 
INACH permits. 

DNA Amplification and Sequencing 

Genomic DNA was extracted from blood samples employing 
saline method [36] . Hypervariable Region I of mtDNA (HVR-I) 
was amplified employing primers from Roeder et al. [37]: 
tRNAGlu (5 ' -CCCGCTTGGCTTYTCTCCAAGGTC-3 ') and 
AH530 (5'- CTGATTTCACGTGAGGAGACCG-3'). Each 
25 u\l PCR reaction contained 0.4 um each primer, 200 u,m each 
dNTP, 1.5 mil MgCl 2 , 1 u.1 10X PCR reaction buffer and 1 U 
Taq Polymerase Invitrogen. Thermocycler conditions were an 
initial denaturing at 95°C for 5 min, followed by 38 cycles of 95°C 
for 30 s, 55°C for 45 s and 72°C for 15 min, and a final extension 
at 72°C for 7 min. 

Intron 7 of the fi-fibrinogen gene (fijtbintT) was amplified using 
the primers Fib7U (5 ' -C AGGAC AATGACAATTC AC-3 ') and 
Fib7L (5'-GTAGTATCTGCCATTAGG-3') [38]. The PCR 
reaction mixture contained 0.4 |j,m each primer, 200 u.m each 
dNTP, 1.5 mM MgCl 2 , 1 ul 10X PCR reaction buffer and 1 U 
Taq Polymerase Invitrogen. Thermal conditions were an initial 
denaturing at 95°C for 10 min, followed by touchdown of 10 
cycles of 95°C for 15 s, 60-50°C for 30 s and 72°C for 45 s, a 
second stage of 35 cycles of 95°C for 15 s, 50°C for 30 s and 72°C 
for 45 s and a final extension phase of 72°C for 30 min. All PCR 
products were visualized in 2% agarose gels and sequenced using 
the service of Macrogen, South Korea. Posterior analysis and 
editing of sequences employed Proseq v 2.0 [39] obtaining 371 bp 
for HVR-I and 1014 bp for [Sfibint7. Sequences were deposited in 
GenBank under accession numbers KF717669-KF7 17743 for 
mtDNA HVR-I, and KF803989-KF803995 for nDNA |3fibint7. 

Data Analysis 

The nuclear marker Pfibint7 haplotype reconstruction was 
developed employing the Phase module in DnaSP v. 5.0 [40]. 
Genotype reconstruction was based on the algorithms of PHASE 
v.2.1 [41,42], using coalescent Bayesian methodology to infer 
haplotypes. Standard diversity indexes (K, S, Hd, n and II) were 
calculated for both HVR-I and pTibint7 using DnaSP v. 5.0 [40]. 
Genealogical relationships between haplotypes were established 
employing the median joining algorithm implemented in Network 
v. 4.61 1 [43]. Assignment weights for frequent mutation steps (>6 
times) were changed from 10 to 5 to eliminate loops and simplify 
the tree. 
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To evaluate population genetic structure we estimated F st and 
0 sf> using Arlequin v.3.5.1.3 [44]. Phylogeographic structure was 
evaluated employing Permut [45,46], considering 10,000 permu- 
tations. Isolation by distance was evaluated with the Mantel test 
(10,000 permutations) implemented in Arlequin v.3.5.1.3 [44]. 
Spatial structure was assessed using Geneland 4.0.2 [47,48] 
occupying a correlated analysis and 1,000,000 iterations. 

Mutation-drift equilibrium was tested employing Tajima's D 
and Fu's Fs indexes [49,50] in Arlequin v.3.5.1.3 [44]. Due to the 
low genetic diversity present in (3fibint7, demographic inference 
analyses were not performed on this marker to avoid overestima- 
tion. Mismatch distribution was reconstructed to evaluate past 
population changes, estimating both demographic and expansion 
parameters under the instantaneous growth model [51]. Expan- 



sion time was estimated from TJ = 2 ut, u. : mutation rate 
substitution/site/year employing a 55% mutation rate for HVR- 
I [52]. Effective population size and population demographic 
history of penguin populations were estimated using Bayesian 
inference employing BEAST v. 1 .7.4 [53] . The evolution model for 
the sequence dataset was selected employing JModelTest under 
Bayesian inference criteria (BIC) [54]. For HVR-I, the HKY+I+G 
model was selected for the data set. Selection of the inference 
method was performed using Bayes Factor [55,56], which selected 
the Bayesian skyline plot for P. papua with a relaxed molecular 
evolution clock and exponential distribution using 0.55 s/s/ma 
mutation rate. Graphic representation was developed in Tracer v 
1.5 [57]. 
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Results 



We identified a total of 75 haplotypes (n = 150) for HVR-I of P. 
papua (371 bp), with high levels of overall haplotype diversity 
(h — 0.982), ranging from 0.924 to 0.975 between localities, and a 
reduced level of nucleotide diversity (71 = 0.01288), ranging from 
0.01009 to 0.01599. For Pfibint7 (1014 bp), a total of 5 haplotypes 
(n= 176) were found, with haplotype diversity of 0.492 and 
reduced nucleotide diversity (it = 0.00091) as expected for a 
nuclear marker with low mutation rate (Table 1). 

The median joining network for mtDNA HVR-I reflects the 
high levels of genetic diversity; there were shared haplotypes 
between localities. However, only one shared haplotype was found 
between the southernmost GGV and other localities (Figure 2a). 
The median joining network for jifibint7, even if there were low 
levels of genetic diversity, shows the existence of shared haplotypes 
between all localities (Figure 2b). 

Fu's Fs and Tajima's D test for HVRT for P. papua was not 
significant (Fs= -7.59286, p = 0.0206; D= - 1.01 107, p = 0. 1644, 
respectively). Mismatch distribution (Figure 3) showed a unimodal 
shape, suggesting a recent expansion (^ = 4.924), which was 
estimated at about 12,000 BP. The Bayesian skyline plot for P. 
papua (Figure 4) also suggests a signature of recent expansion 
around 1 3,000 BP with the Most Recent Common Ancestor time 
(MRCA) around 21,500 BP. The increase in effective population 
size (Ne) for P. papua is about one order of magnitude. 

Spatial genetic analysis to infer subpopulations (K) indicated an 
absence of clusters (K= 1). However, a reduced but significant 
population genetic structure was observed between pairwise 
populations for HVRT for P. papua, with F st values ranging from 
0.027 to 0.052 and 0„ ranging from 0.041 to 0.102 (Table 2). The 
only non-significant values of F st and 0 S , were found between the 
Ardley Island and Admiralty Bay sites, located only 30 km apart in 
King George Island. We did not find a significant rignature of 
isolation by distance (r= -0.046770, p = 0.5015). 

Discussion 

Molecular data support the hypothesis that P. papua populations 
were affected by Pleistocene climate changes. This can be 
evidenced by high levels of genetic diversity, an absence of 
population expansion signature from Fu's Fs and Tajima's D, 
whereas the expansion signature detected by Bayesian skyline plot 
can be associated to a recolonization scenario (see below). We 
found a signature of population expansion for P. papua after the 
LGM (12,000-13,000 BP), and the MRCA around 21,500 BP. 
The date of expansion calculated is coincident with the retreat of 
the ice sheet cover in the maritime Antarctic (Antarctic Peninsula 
and adjoining islands, [58]). Nevertheless, fossil evidence suggests 
that P. papua would have established colonies in South Shetland 
Islands much later, around 1,100 BP during the late Holocene 
[34] . However fossil records for P. papua axe very scarce and it is 
possible that older evidence has not yet been recovered. There is 
also the possibility that our expansion time could be overestimated 
due to the great amount of genetic diversity, which can be 
attributed to a colonization of maritime Antarctic from a large 
number of migrant individuals from different refuges. Dinechin et 
al. [59] detected the existence of three different clades in P. papua: 
one associated with sub-Antarctic islands near the Indian Ocean 
(Kerguelen), another associated with the Falkland Islands and the 
third associated with maritime Antarctic and peri-Antarctic 
Islands (i.e. South Georgia, Gough, Sandwich Islands [58]). The 
existence of a unique linage in the maritime Antarctic suggests that 
during the LGM more than one colony remained near the 
Antarctic Peninsula, in the East Antarctic Peninsula and some 
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Figure 2. Median Joining Network of mtDNA HVR-I (a) and pfibint7 for P. papua (b). GGV: Gabriel Gonzalez Videla Base. 
doi:1 0.1 371 /journal.pone.0095375.g002 



maritime and peri-Antarctic island refugia, as has been proposed 
for another Antarctic species [28]. Contrary to the expected for 
this scenario a severe signature of bottleneck was not detected, and 
all populations exhibited high levels of genetic diversity. These 
results may indicate that the colonization of deglaciated areas 
involved a large number of individuals from multiple locations or 
refuges, rather than through successive founder effects as has been 
proposed for many Antarctic species [33,60]. Therefore, our 
results suggest that the species may have survived in peri-Antarctic 
refugia such as South Georgia and North Sandwich Islands and 
recolonized the Antarctic Peninsula and South Shetland Islands 
after the ice sheet retreat. Recent colonization from multiple 



sources is also supported by the low levels of genetic structure 
found in our sampling area as well as the absence of isolation by 
distance among sites. 

Although low, mtDNA HVR-I genetic differentiation among 
localities was significant between most of the pairwise compari- 
sons, except between the two closest localities in King George 
Island (Admiralty Bay and Ardley Island); it was not possible to 
detect any genetic structure with nDNA |3fibint7. However, this 
discrepancy is likely related to the difference in substitution rates 
and diversity between the markers [61,62]. Alternatively, differ- 
ences in genetic structure between nuclear (biparental markers) 
and mitochondrial markers (maternal lineage) can result of 
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Figure 3. Mismatch Distribution for mtDNA HVR-I for P. papua. 
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differential migration or dispersal rates between males and 
females. Philopatric behavior has been well established in 
pygoscelid penguins; adults breed in or near their natal colony 
[63] . For instance, P. papua displays a relatively high degree of site 
fidelity, between 89-100% in South Georgia [64] and 60% in 
King George Island [65]. Unfortunately, no information is 
available about differences in philopatry between sexes in this 
species. It is thus not possible to support a female-biased philopatry 
hypothesis to explain the discrepancy between mitochondrial and 
nuclear markers. Moreover, in the closely related P. adeliae, [66] 
found that nest site fidelity is higher for males (98.9%) than for 
females (65.3%). Another study found also that the majority of 
males (62%) returned to the same nest during the 4 years of study, 
while only 29% females returned to the same nest [67]. Although 



philopatric behavior is not concordant with a low genetic structure 
found for mtDNA, which is a maternal lineage molecular marker, 
the colonizing ability of P. papua could be an explanation. Gentoo 
penguins are apparendy excellent colonizers of new breeding 
territory, especially at their southern range boundary [68]. 
Moreover, range expansion erases phylogeographic structure, as 
suggested for P. adeliae [31]; this pattern was observed for other 
Antarctic bird species. The south polar skua has undergone 
demographic and/or spatial expansions and is considerably less 
phylogeographically structured than the brown skua [69]. 
Additionally, it has been proposed that philopatric conduct may 
be related to stable environmental conditions [70]. Moreover, 
Dugger et al. [71] found that environmental perturbations can 
increment dispersal rate among Adelie penguins; they reported 
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Figure 4. Skyline Plot mtDNA HVR-I for P. papua. 
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Table 2. Population genetic structure between locations of P. papua. 
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that in periods of high environmental adversity penguins tend to 
increment their movement rates in all colonies and in all 
directions. Consistent with this idea of high dispersion in the face 
of climate constrains, Korczak-Abshire et al. [72] studied AFLP 
for P. antarcticus in two localities from King George Island, finding 
no or only weak population genetic structure. Additionally, 
Roeder et al. [73] using microsatellites also found reduction in 
and even absence of genetic diversity for P. adeliae around the 
Antarctic. Considering this information, it is expected that low 
values of genetic structure may be associated with a relaxing of 
philopatric behavior in P. papua, allowing movement and 
establishment of colonies of individuals from different source 
colonies.Further studies along the entire P. papua distribution using 
multilocus markers will help to completely understand the 
evolutionary history of this species, and the results of biparental 
variable molecular markers (e.g. microsatellite loci) can be 
compared with our results of population genetic structure from 
mtDNA (maternal lineage) to further understand the species 
behavior. 

Conclusions 

Understanding how past climate change affected species 
demography could help to elucidate the impact of current global 
warming and to predict future trends in species. The study of P. 



papua demography allowed us to detect that in the past P. papua 
maintained a large genetic diversity associated with large effective 
population size. Although it has been described as a sub-Antarctic 
species, our results suggest that the Gentoo penguin probably 
persisted in peri-Antarctic refugia during the LGM, being able to 
(re)colonize the maritime Antarctic during warmer periods (such 
"penguin optimum" period, Baroni and Orombello, 1994). 
According to this, P. papua may be considered as a resilient species 
which may have survived the extreme cold temperature during the 
LGM, and is obtaining benefits from the recent climate change in 
the maritime Antarctic. Over the last 30 years Gentoo penguin 
populations have been increasing in the WAP [18,74], colonizing 
and establishing new colonies with fast growth rates along the 
southern boundary [20], and it is rapidly adapting to the shifts in 
the marine food web due to its opportunistic foraging [68]. It is 
expected that in the future P. papua will be a dominant species, at 
least in the Antarctic Peninsula, as temperatures increase in the 
Antarctic. 
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